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Fluctuations in nucleon positions can affect the spatial eccentricity of the overlap zone in nucleus- 
nucleus collisions. We show that elliptic flow should be scaled by different eccentricities depending 
on which method is used for the flow analysis. These eccentricities are estimated semi-analytically. 
When V2 is analyzed from 4-particle cumulants, or using the event plane from directed flow in a 
zero-degree calorimeter, the result is shown to be insensitive to eccentricity fluctuations. 

PACS numbers: 25.75.Ld, 24.10.Nz 



1. Introduction 

Elliptic flow, V2, is one of the key observables in 
nucleus-nucleus collisions at RHIC 1]. It originates 
from the almond shape of the overlap zone (see Fig. 
which produces, through unequal pressure gradients, an 
anisotropy in the transverse momentum distribution 0, 
the so-called V2 = (cos 20} , where 0's are the azimuthal 
angles of the detected particles with respect to the reac- 
tion plane. 

Preliminary analyses of V2 in Cu-Cu collisions at RHIC 
HHH, presented at the QM'2005 confer cncc, reported 
values surprisingly large compared to theoretical expec- 
tations, almost as large as in Au-Au collisions. It was 
shown by the PHOBOS collaboration Q that fluctua- 
tions in nucleon positions provide a natural explanation 
for this large magnitude. The idea is the following: The 
time scale of the nucleus-nucleus collision at RHIC is so 
short that each nucleus sees the nucleus coming in the op- 
posite direction in a frozen configuration, with nucleons 
located at positions whose probabilities are determined 
according to the nuclear wave function. Fluctuations in 
the nucleon positions result in fluctuations in the almond 
shape and orientation (see Fig. and hence in larger 
values of V2 . 

In this Letter, we discuss various definitions of the ec- 
centricity of the overlap zone. We show that estimates 
of V2 using different methods should be scaled by appro- 
priate choices of the eccentricity. We then compute the 
effect of fluctuations on the eccentricity semi-analytically 
to leading order in 1/iV, where N is the mean number 
of participants at a given centrality. A similar study was 
recently performed by S. Voloshin on the basis of Monte- 
Carlo Glauber calculations 6]. 

2. Eccentricity scaling and fluctuations 

Elliptic flow is determined by the initial density pro- 
file. Although its precise value depends on the detailed 
shape of the profile, most of the relevant information is 
encoded in three quantities: 1) the initial eccentricity of 
the overlap zone, e, which will be defined more precisely 
below; 2) the density n, which determines pressure gradi- 
ents through the equation of state (by density, we mean 
the particle density, n, at the time when elliptic flow de- 



FIG. 1: Schematic view of a collision of two identical nuclei, 
in the plane transverse to the beam direction (z-axis). The 
X- and j/-axes are drawn as per the standard convention. The 
dots indicate the positions of participant nucleons. Due to 
fluctuations, the overlap zone could be shifted and tilted with 
respect to the (a;, y) frame, x' and y' are the principal axes 
of inertia of the dots. 



velops; this time is of the order of the transverse size R. 
Quite remarkably, the density thus defined varies little 
with centrality, and has almost the same value in Au-Au 
and Cu-Cu collisions at the same colliding energy per 
nucleon 7J); 3) the system transverse size R, which de- 
termines the number of collisions per particle. V2 scales 
like e for small e, that is, V2 — £f{n, R). 

This proportionality relation is only approximate. 
However, hydrodynamical calculations (3 show that it is 
a very good approximation in practice for nucleus- nucleus 
collisions. Eccentricity scaling holds for integrated flow 
as well as for the differential flow of identified particles. 
In the latter case, the function f{n,R) also depends on 
the mass, transverse momentum and rapidity of the par- 
ticle. 

Eccentricity scaling of V2 is generally believed to be a 
specific prediction of relativistic hydrodynamics. In the 
form above, the scaling is expected to be more general: it 
does not require thermalization, as implicitly assumed by 
hydrodynamics. If thermalization is achieved, that is, if 
the system size R is much larger than the mean free path 
A, then the scaling is stronger: i^/e no longer depends 



on R, but only on the density n 
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The standard definition of the eccentricity is 

~ 



(1) 



where (x, y) is the position of a participant nucleon in 
the coordinate system defined in Fig. ^ Throughout this 
paper, (...) denotes an ensemble average: here, it means 
an average over participant nucleons and over many coUi- 
sion events of the same impact parameter. This standard 
eccentricity apphes to most hydrodynamic calculations: 
indeed, most hydrodynamic calculations (with the excep- 
tion of Ref. la) use smooth, event-averaged initial con- 
ditions, with an initial entropy density proportional to 
the density of participants ilfl] . The charged multiplicity 
then scales like the number of participants, as observed 
experimentally. (Adding a component proportional to 
the number of binary collisions, as argued in [ILJ . does 
not change the eccentricity significantly) . 

It was recently argued [l^] that the Color Glass Con- 
densate picture of heavy-ion collision leads to a different 
definition of the eccentricity, which may be significantly 
larger than the standard eccentricity. This interesting 
possibility will not be considered further here. 

Now, because of the event-b y-e vent fluctuations in the 
participant nucleon positions [l^, the eccentricity driv- 
ing elliptic flow in a given event is that defined by the 
principal axes {x',y') of the distribution of participant 
nucleons, see Fig. This "participant eccentricity" Epart 
can be written as [3 



-part 



where 



{xy} - {x}{y}, 



(2) 



(3) 



and {...} denotes the average over all participants in one 
collision event (sample average). Our basic assumption 
in this paper is that in each event, the elliptic flow V2 is 
proportional to the participant eccentricity epart- 

Experimentally, elliptic flow is analyzed by selecting 
events in a centrality class. The quantities which drive 
elliptic flow, namely, the density n, the transverse size 
i?, the eccentricity Epart, fluctuate from one event to the 
other. This causes dynamical fluctuations of W2. The ef- 
fect of imp act parameter fluctuations was carefully stud- 
ied in [iJI and was shown to be small. In this paper, 
we focus on fluctuations in the positions of participant 
nucleons. It will be shown below that fluctuations in the 
eccentricity Epart dominate over fluctuations in size and 
density. 

3. To each method its own eccentricity 

Since the reaction plane is not known exactly on an 
event-by-event basis, V2 is measured indirectly using az- 
imuthal correlations. Several methods have been used. 



which yield different estimates of V2 ■ We argue that these 
estimates are affected in different ways by fluctuations of 

• The event-plane method [Tsl l has been implemented 
by the STAR 1], PHOBOS Te^ and PHENIX [l^ 
collaborations at RHIC. The event plane is an es- 
timate of the reaction plane; it is deflned as the 
plane spanned by the collision axis and the ma- 
jor axis of the ellipse formed by the transverse 
momenta of outgoing particles. It corresponds to 
the x' axis defined by the participants (Fig. 
up to statistical fluctuations which are taken care 
of by the analysis (this is the so-called "event- 
plane resolution"). Outgoing particles are then 
individually correlated to this event plane. The 
corresponding estimate of V2 will be denoted by 
W2{EP2}. Another method determines V2 from two- 
particle azimuthal correlations between outgoing 
particles 0. The corresponding estimate is usu- 
ally denoted by W2{2}, and is deflned by an equa- 
tion of the type W2{2} — \J (cos 2(01 — 02))- Both 
methods are essentially equivalent (see Sec. III.C.2 
of Ref. If ■'^2 fluctuates, both yield the rms 

value of V2 El: i;2{2} ~ W2{EP2} = yfWl- Since 
V2 in each event scales with Epartj one expects 
W2{2|~ W2{EP2} oc e{2}, where £{2} is defined 

by [ia 



r{2} ^ {^). 



(4) 



This scaling differs from that proposed in Ref. jj], 

V2 OC (Spart)- 

• Four- par ticle cumulants of azimuthal correla- 
tions [23 can be used to reduce the bias from non- 
flow correlations. The corresponding estimate of V2 
is denoted by t'2{4}. It involves a combination of 
2-particle and 4-particle correlations, i.e., the 2nd 
and 4th moments of the distribution of W2, and 
(wl). Assuming again that V2 scales with £part in 
each event, one obtains W2{4} oc e{4}, where e{4} 
is defined by [l^ : 



£{4}^(2(4,,,)'-(£4p,,,)) 



1/4 



(5) 



• An alternative method determines the event-plane 
from directed flow [2ll | , determined in a zero-degree 
calorimeter (ZDC) which detects spectator neu- 
trons. Particles in the central rapidity region are 
then correlated to this event plane. This last es- 
timate of V2 is denoted by W2{ZDC} j5|. It differs 
from the previous ones in that the reference direc- 
tion (the event plane) is determined by spectator 
neutrons from the projectile, rather than by partic- 
ipants. The direction deflned by spectator neutrons 
is the X-axis, not the cc'-axis (see Fig.Pl, up to fluc- 
tuations in spectator positions which are taken care 
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of by the analysis. Therefore, the relevant eccen- 
tricity for this analysis is the reaction-plane eccen- 
tricity 



(6) 



and U2{ZDC} should scale correspondingly like 

V2{ZBC} cx (£Rp) . (7) 

This estimate of V2 is also unbiased by nonflow correla- 
tions because it involves a 3-particle correlation (instead 
of a 2-particle correlation in the standard event-plane 
method), and also because the ZDC calorimeter has a 
wide rapidity gap with the central rapidity detector. 

4. Computing the fluctuations 

We now compute the eccentricities entering Eqs. 
© and lO, namely, e{2}, e{4} and (erp), to leading 
order in the fluctuations. 

To that end, we write each sample average (over one 
event) entering Eq. ^ as the sum of an ensemble average 
and a fluctuation, e.g., {x^} = {x^) + ^x^, where 8^1.2 
is the fluctuation. We choose the coordinate system in 
Fig. n where {x) ~ {y) = (xy) ~ 0. Substituting in 
Eq. (PI and retaining terms to second order in the J's 
gives 



-part 



Ss + 



Sy2 ~ 5x 



Sy2 + 6y - 5^ 



6\ - 5^ 

xy y-^ X 



(r2 



(8) 



where (r"^^ = l^x"^ +2/^)- The reaction plane eccentric- 
ity Erp is given by a similar expression, except for the 
fifth term on the right-hand side (rhs) proportional to 5'^y 
which does not appear in SRp. The 2nd and 3rd terms on 
the rhs of Eq. ||SJ| arc linear and the remaining 5 terms are 
quadratic in fluctuations. In an obvious notation Eq. (jHl 
can be rewritten as 

2 5 

£part - + ^ A,<5f ) + M^'' : (9) 
i=l i=l 

where 5^'^' are linear and Sf'^ are quadratic in J's. erp 
is given by a similar expression, with one less quadratic 
term. It is now straightforward to derive the expressions 

for (epart) , (epart) : (4art)> etc. Wc get 

(epart) = es+(-B,5f^) 

{eU) = et + %el{{A,5f^f)+Ael{B..5f')^ 

(10) 



where we have retained terms to second order in the 5's 
and used ^^i^^^ = 0; summation over the repeated index 
is understood. 

The averages on the rhs of Eqs. H1U|) are easily com- 
puted by using the following identity, which holds for N 
independent participant nuclcons 



if 9) if) (g) 
N 



(11) 



where / and g are any functions of x and y. We get 

e{2r = el+ 

(r^) ( , (r4cos2,^)\ , ^ 

I 1 + 3e? + Ae.^^rr-^ , (12) 



iV(r2) 



(r^ cos 40) 



and 



(£Rp) = £s + 



(r^ cos 2q 



(13) 



(14) 



where (r, (p) are the polar coordinates in the (x, y) plane. 

The number of participants, TV, and the moments of 
their distribution, (?''*), (r'' cos 20), . . . , are easily com- 
puted numerically in the Glauber model with Woods- 
Saxon distribution of participants. The parameters of 
our Glauber calculation are the same as in Ref. [lTI| . 
with one minor difference: we neglect fluctuations in the 
number of participants, i.e., we assume a one-to-one re- 
lation between impact parameter and number of partici- 
pants. Orders of magnitude of the various coefficients in- 
volved are (r^) / {r'^f ^ 1.4-2, (r^ cos 20) / (r^) ~ -e^, 
(r"* cos 40) / (r"') e^. Numerical results for the various 
eccentricities as a function of the number of participants 
are given in Fig. |21for Au-Au and Cu-Cu collisions. 

The standard eccentricity Eg vanishes for the most cen- 
tral as well as for the most peripheral collisions as ex- 
pected. On the other hand, e{2} rises steeply as iV de- 
creases; this is due to the 1/A^ contribution of the fluc- 
tuations, see Eq. (|12l) . This contribution originates from 
the 5y2 — (5^2 and (5^^, terms in Eq. The number of 
participants N being smaller for the Cu-Cu collision than 
for the Au-Au collision, the magnitude of the fluctuations 
is relatively larger in the former case. 

The 4-cumulant e{4} is only slightly larger than £s, 
down to low values of N . This is because the 1/A^ term 
in Eq. H13|) is multiplied by a term of order e^, As a 
result, the relative difference between e{4} and £s is only 
of order 1/A^. Similarly, the reaction plane eccentricity, 
(erp), is almost equal to £s- 

We have considered so far the fluctuations in the ec- 
centricity of the overlap zone. We now briefly discuss 
fiuctuations in the size or the transverse area of the over- 
lap zone, S ~ 2^(7 x<Jy They can be calculated exactly 
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FIG. 2: Standard eccentricity Sa, participant eccentricities 
£{2} and £{4}, and reaction-plane eccentricity (£rp), vs the 
number of participant nucleons for a Au-Au collision (top) 
and a Cu-Cu collision (bottom). 



ill the same way as fluctuations in the eccentricity. The 
result is 




{Sf 

(15) 

Numerically, these fluctuations are found to be practi- 
cally negligible. The reason why eccentricity fluctuations 
are important is that the term in Eq. (|12|l must be 
compared with e^, which is itself a small number. 

5. Discussion 

We have shown that the values of V2 analyzed with 
different methods should be scaled by different eccentric- 
ities: U2{EP2} (standard V2), W2{4} and W2{ZDC} should 
be scaled respectively by e{2}, e{4} and (erp), defined 
in Eqs. (jSJ and l(7|). 

Our most important result is that e{4} and (£rp) are 
almost equal to the standard eccentricity, while £{2} is 
strongly affected by fluctuations for small systems and/or 
peripheral collisions. An important contribution to the 
fluctuations comes from the angle between the x'-axis 



and the x-axis in Fig.^ i.e., from the angle of tilt of the 
participant ellipse relative to the reaction plane. This ef- 
fect was neglected in Ref. 0| and first taken into account 
m Ref. % The ZDC analysis eliminates the effect of the 
tilt angle by measuring the eccentricity in the (x,y) axes; 
in the case of the 4-cumulant analysis, most of the fluc- 
tuations happen to cancel in the subtraction of Eq. 10). 

Our results show that higher-order estimates of ellip- 
tic flow, U2{4} and W2{ZDC}, are not only insensitive to 
nonflow effects, but also, to a large extent, to fluctua- 
tions in the participant eccentricity. This is confirmed 
by transport calculations , and explains the observed 
agreement between U2{4} and U2{ZDC} in Au-Au colli- 
sions 0. 

Fluctuations in the participant eccentricity, on the 
other hand, tend to increase the value of the standard, 
event-plane V2. We have estimated this increase quan- 
titatively, assuming independent nucleons. Within this 
simple model, fluctuations account for less than one half 
of the observed difference between W2{2} and W2{4} in 
Au-Au collisions Q . The remaining difference can be (at 
least partly) ascribed to nonflow effects. Nonfiow effects 
are clearly seen in the different pr-dependences of W2{2} 
and W2{4}, which cannot be explained by fluctuations. 

However, one should be aware that estimates of the 
participant eccentricity are model dependent: in partic- 
ular, Monte-Carlo Glauber calculations 0, generally 
yield higher fluctuations. The difference is the follow- 
ing: in a Monte-Carlo Glauber, each nucleon is modelled 
as black disk of transverse area a, and a nucleon from 
nucleus A is a participant if it overlaps with at least 
a nucleon from nucleus B, and vice- versa. Participant 
nucleons are therefore correlated, and the black-disk ap- 
proximation maximizes these correlations. Due to such 
correlations, which are neglected in our calculation, our 
estimate of the participant eccentricity can be considered 
a lower bound. Quantitatively, correlations can increase 
the effect of fluctuations by a factor up to 2: if the partic- 
ipant nucleons can only be found in pairs of overlapping 
disks, this amounts to replacing N by N/2 (the number 
of pairs) in Eqs. H12I14|I . 

Finally, let us compare our results with the recent 
Monte-Carlo Glauber calculations of Ref. The re- 
sults are in qualitative agreement with ours: e{4} is much 
closer to Es than to £{2} for moderate centralities. For 
large impact parameters, however, £{4} is almost equal 
to £{2}. This means that the participant eccentricity, al- 
though much larger than the standard eccentricity, fluc- 
tuates little from one event to the other. This intriguing 
behaviour could be a consequence of the strong corre- 
lations mentioned above, and deserves further investiga- 
tion. 

To summarize, the elliptic flow scaled by the eccentric- 
ity of the overlap zone, V2/£, is an important observable 
at RHIC as well as LHC, because if it is found to be in- 
dependent of the system size, one has a strong pointer 
toward thermalization. We have discussed various defi- 
nitions of V2 and e and studied how they are affected by 
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fluctuations in nucleon positions. We have shown that 
when V2 is analyzed using 4-particle cumulants or the 
event-plane from directed flow in a ZDC calorimeter, the 
resulting estimate essentially scales with the standard ec- 
centricity, and is insensitive to the fluctuations in the 
participant eccentricity considered by PHOBOS 0. 



Wang, G. Roland, M. Nardi, S. Manly and A. Poskanzer 
for discussions. 
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